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Abstract 

As was shown in an earlier paper pQ, solids dispersed in a drying drop migrate to 
the (pinned) contact line. This migration is caused by outward flows driven by the loss 
of the solvent due to evaporation and by geometrical constraint that the drop maintains 
an equilibrium surface shape with a fixed boundary. Here, in continuation of our earlier 
paper we theoretically investigate the evaporation rate, the flow field and the rate 
of growth of the deposit patterns in a drop over an angular sector on a plane substrate. 
Asymptotic power laws near the vertex (as distance to the vertex goes to zero) are 
obtained. A hydrodynamic model of fluid flow near the singularity of the vertex is 
developed and the velocity field is obtained. The rate of the deposit growth near the 
contact line is found in two time regimes. The deposited mass falls off as a weak power 
7 of distance close to the vertex and as a stronger power (5 of distance further from the 
vertex. The power 7 depends only slightly on the opening angle a and stays between 
roughly — 1/3 and 0. The power (5 varies from —1 to as the opening angle increases 
from to 180 degrees. At a given distance from the vertex, the deposited mass grows 
faster and faster with time, with the greatest increase in the growth rate occurring at 
the early stages of the drying process. 

PACS: 47.55. Dz — Drops and bubbles; 68.03.Fg — Evaporation and condensation; 
81.15.-z — Methods of deposition of films and coatings; film growth and epitaxy. 

1 Introduction 

The problem of the so-called "coffee-drop deposits" has recently aroused great interest [H 
El IU Ej- The residue left when coffee dries on the countertop, mineral rings left on washed 
glassware, banded deposits of salt on the sidewalk during winter and enhanced edges in water 
color paintings are examples of the variety of physical systems understood by "coffee-drop 
deposits" terminology. Understanding the process of drying of such solutions is important 
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for paint manufacturing, protein crystallography, production of nanowires, patterning of a 
surface, DNA stretching in a flow, and other scientific and industrial applications. 

The theory of the solute transfer in such physical systems has been developed in works 
of Deegan et al. |U El El- In this theory, the contact line of a drop of liquid is pinned 
during drying process. While the highest evaporation occurs at the edges, the bulk of the 
solvent is concentrated closer to the center of the drop. In order to replenish the liquid 
removed by evaporation at the edge, a flow from the inner to the outer regions must exist 
inside the drop. This flow is capable of transferring all of the solute to the contact line 
and thus accounts for the strong contact-line concentration of the residue left after complete 
drying. The theory of Deegan et al. [1J is very robust since it only requires the pinning of 
the edge during drying, which can occur in a number of possible ways (surface roughness, 
chemical heterogeneities etc), and it is independent of the nature of the solute. It accounts 
quantitatively for the experimentally observed phenomena at least in the case of geometry 
analyzed in [T]. However, only the simplest case of a round drop was analytically solved by 
Deegan et al. 

Here we consider a complementary problem of a solute drop drying over an angular region 
(Fig. [IJ. An arbitrary boundary line can be represented as a sequence of smooth segments, 
which can be approximated by circular arcs, and fractures, which can be approximated by 
angular regions. Thus, knowledge of analytical solution for both circular (solved by Deegan 
et al.) and angular (considered here) boundary shapes fills out the quantitative picture of 
solute transfer and deposit growth for an arbitrary drop boundary. Keeping this purpose 
in mind, we specify only one boundary of the drop (the vertex and the sides of the angle) 
leaving the remainder of the boundary curve unspecified. Such approach turns out to be 
sufficient to determine the universal features of the solution, and it allows us to find all the 
important singularities as power laws of distance from the vertex of the angle. 

Our motivation for focusing on droplets over an angular sector also arises from the distinc- 
tive deposit patterns observed in Deegan's experiments (Fig.|2J). These patterns are produced 
as the contact line retracts (down on the figure) discontinuously, switching between pinned 
and moving states. We expect that the knowledge of the solution for a single pinned angle 
can boost understanding of these distinctive dynamically-produced patterns. 

The flow pattern explored in Deegan's work |3| is a new type of capillary flow that 
depends sensitively on the perimeter shape. The range of flows and deposition profiles ob- 
tainable by this mechanism remains virtually unexplored. Our study aims to explore the 
extremes of what behavior can be obtained by varying the perimeter shape. 

As a result of our study, we find that the solution in the angular-sector geometry is in a 
number of ways different from the solution in the circular geometry. In particular, three time 
regimes for the deposit growth are found in the angular case compared to the two regimes in 
the circular one. A clear-cut signature of the existence of these new regimes is also provided. 
In addition, the new geometry possesses an extra free parameter — the opening angle of the 
sector, — and dependence of all the universal exponents in the power laws for all quantities 
on this extra parameter yields the possibility to control the deposition patterns by simply 
choosing a proper opening angle. Thus, we obtain results that are both universal (i.e. do not 
depend on any physical properties of the constituents) and at the same time dependent on 
the geometry in a predictable fashion. These and other useful properties of our results are 
discussed at the end of this manuscript. 
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(b) 

Figure 1: (a) A water droplet with a sector-shaped boundary on the plane substrate (side 
view), (b) The same droplet pictured from another point (top view). Black lines are the 
grooves on the substrate necessary to "pin" the contact line. (Photos by Itai Cohen.) 
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Figure 2: Deposit patterns observed in experiments by Deegan at al. These patterns were 
left after the contact line moved through the area shown. The contact line was retreating 
down on this image and the solute was below the contact line. (Photo by Robert Deegan.) 



In the following section we first consider the physical assumptions and the mathematical 
formulation of our theory, then calculate the equilibrium surface shape of the drop and the 
rate of evaporation from the surface of the drop, and finally obtain the velocity field inside 
the drop. At the end of the section we describe the trajectories of the particles and provide 
the results for the solute transfer to the contact line and the deposit growth in different time 
regimes. In the last section we discuss achievements and limitations of our theory and explore 
the avenues for further study. 

2 Theory 

System and assumptions. We consider a droplet of solution on a horizontal surface 
bounded by an angle a in the plane of the substrate (Fig. [TJ. We assume that the droplet 
is sufficiently small so that the surface tension is dominant, and the gravitational effects can 
be safely neglected. At the same time, we do not assume that the contact angle between the 
liquid-gas interface and the plane is constant along the boundary line on the substrate, nor 
do we assume it is constant in time. To achieve an angular boundary, the substrate must have 
scratches, grooves or other inhomogeneities (sufficiently small compared to the dimensions of 
the droplet), which pin the contact line. A strongly pinned contact line can sustain a wide 
range of contact angles; the angle is not fixed by the interfacial tensions as it is on a uniform 
surface. 

The use of the cylindrical coordinates (r, 0, z) is most natural in this problem, so that 
the angle occupied by the liquid on the substrate is < r < oo and —a/2 < (j) < a/2, 
and the coordinate normal to the substrate is z (Fig. EJ). The geometry of this problem is 



4 



bisector 




®Z 

Figure 3: Geometry of the problem. The plane of the figure coincides with the substrate. 

much more complicated than the one of the round-drop case solved earlier. We consider 
a two-dimensional object — the angular sector — in a three-dimensional space, and the 
main complication arises from the fact that the symmetry of the object does not match the 
symmetry of any simple orthogonal coordinate system in that space. In particular, even 
the solution of the Laplace equation (needed below) requires introduction of the special 
coordinate system (the so-called conical coordinates, or the orthogonal coordinates of the 
elliptic cone) with heavy use of the Jacobi elliptic functions. Similarly, a separate research [2] 
was required to find the equilibrium surface shape of the drop in this geometry. Thus, given 
this complex geometry, we limit our task to determining only the power-law scaling for most 
quantities, and this task proves to be sufficiently challenging by itself. 

We describe the surface shape of the drop h(r, <ft) by local mean curvature H that is 
spatially uniform at any given moment of time, but changes with time as droplet dries. 
Ideally, the surface shape should be considered dynamically together with the flow field 
inside the drop. However, as we show in the Appendix, for flow velocities much lower than 
the characteristic velocity v* = a/3r] (a is surface tension and rj is dynamic viscosity), which 
is about 24 m/s for water under normal conditions, one can consider the surface shape 
independently of the flow and use the equilibrium result at any given moment of time for 
finding the flow at that time. The equilibrium surface shape h(r, 0) for the drop over an 
angular region was found in our earlier paper [2]. 

In order to determine the flow caused by evaporation, one needs to know the flux profile of 
liquid leaving each point of the surface by evaporation. The functional form of the evaporation 
rate J(r, 0) (defined as evaporative mass loss per unit surface area per unit time) depends on 
the rate-limiting step, which can, in principle, be either the transfer rate across the liquid- 
vapor interface or the diffusive relaxation of the saturated vapor layer immediately above 
the drop. As in the work of Deegan et al, we assume that the rate-limiting step is diffusion 
of liquid vapor (Fig. HJ) and that evaporation rapidly attains a steady state. Indeed, the 
transfer rate across the liquid- vapor interface is characterized by the time scale of the order 
of 10 -10 s, while the diffusion process has characteristic times of the order of R 2 /D (where 
D is the diffusion constant for vapor in air and R is a characteristic size of the drop), which 
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Figure 4: The rate-limiting process for evaporative mass loss. It is the diffusion of saturated 
vapor just above the interface rather than the transfer across the interface. 



is of the order of seconds for water drops under typical drying conditions. 

We consider slow flows, i.e. flows with low Reynolds numbers (also known as "creeping 
flows"). This amounts to the neglect of the inertial terms in the Navier-Stokes equation. 
We also employ the so-called "lubrication approximation" . It is essentially based on the two 
conditions reflecting the thinness of the drop and resulting from the separation of the vertical 
and horizontal scales. One is that the pressure inside the drop p does not depend on the 
coordinate z normal to the substrate: d z p = 0. The other is related to the small slope of 
the free surface |V/t| <C 1, which is equivalent to the dominance of the ^-derivatives of any 
component Ui of flow velocity u: d z Ui ^> d s Ui (index s refers to the derivatives with respect 
to any coordinate in the plane of the substrate). The lubrication approximation is a standard 
simplifying procedure for this class of hydrodynamic problems [HI IE] • 

Having formulated physical assumptions intrinsic to the theory, we are now in position 
to formulate its main ideas. 

Equations. The essential idea behind the theory is that a pinned contact line entails 
fluid flow toward that contact line, since the rate of evaporation is the highest at the edge of 
the drop while the most of the liquid is concentrated away from it. The "elasticity" of the 
liquid-air interface fixed at the contact line provides the force driving this flow. 

To develop this idea mathematically, we ignore for a moment any solute in the liquid. 
Once the flow is found, one can track the motion of the suspended particles, since they are 
just carried along by the flow. We define depth-averaged flow velocity by 



where u s is the in-plane component of the local three-dimensional velocity u. Then we write 
the conservation of fluid mass in the form 



where t is the time, p is the density of the fluid and quantities h, J and v are functions 
of r, and t. (We will drop the Vh part of the second term everywhere in the following 
since it is always small compared to unity, as it will be seen from the expression for h(r, <fi) 
below.) This equation represents the fact that the rate of change of the amount of fluid in a 
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Figure 5: Conservation of mass: the liquid-vapor interface lowers exactly by the amount of 
fluid evaporated from the surface plus the difference between the outflow and the influx of 
fluid from the adjacent regions. 



volume element (column) above an infinitesimal area on the substrate (third term) is equal 
to the negative of the sum of the net flux of liquid out of the column (first term) and the 
amount of mass evaporated from the surface element on top of that column (second term); 
Fig. illustrates this idea. Thus, this expression relates the depth-averaged velocity field 
v(r,(f),t) to the liquid-vapor interface position h(r,<f),t) and the evaporation rate J(r,<ft,t). 
However, this is only one equation for two variables since vector v has two components in the 
plane of the substrate. Moreover, as it was pointed out earlier, while the evaporation rate 
J is indeed independent of flow v, the free-surface shape h should in general be determined 
simultaneously with v. Thus, there are actually three unknowns to be determined together 
(h, v r and v^), and hence two more equations are needed. In order to find these equations 
we will employ the lubrication approximation. 

We start with the Navier-Stokes equation with inertial terms omitted (low Reynolds 
numbers): 

Vp = r/V 2 u, (3) 

where p is the fluid pressure, rj is the dynamic viscosity and u is the velocity. Applying 
lubrication-approximation conditions d z p = and d z Ui 3> d s Ui, we arrive at the simplified 
form of this equation 

V s p = r]d zz u S} (4) 

where index s again refers to the vector components along the substrate. From now on we 
will suppress the subscript s at the symbol of nabla-operator, and will assume for the rest of 
this paper that this operator refers to the two-dimensional vector operations in the plane of 
the substrate. Solution to the above equation with boundary conditions 

u s\ z=0 = and d z u s \ z=h = (5) 

yields 

X7r> I r 2 \ 

(6) 
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or, after vertical averaging (Q), 

h 2 

v=--Vp. (7) 

This result is a variant of the Darcy's law jHUHl- Note that since a curl of a gradient is always 
zero and 77 is a constant, the previous equation can be re- written as 

Vx(£)=0. (8) 

This condition is analogous to the condition of the potential flow (V x v = 0), but with a 
quite unusual combination of the velocity and the surface height v/h 2 in place of the usual 
velocity v. 

Relation (|7J) provides the two sought equations in addition to the conservation of mass (j2J). 
However, it contains one new variable, pressure p, and hence another equation is needed. This 
last equation is provided by the condition of the mechanical equilibrium of the liquid-air 
interface (also known as the Young-Laplace equation) relating the pressure and the surface 
shape: 

P = -2H<7 +Patm- (9) 

Here p a tm is the atmospheric pressure, a is the surface tension and H is the mean curvature 
of the surface, uniquely related to the surface shape h by differential geometry. Note that 
this expression is independent of both the conservation of mass (0) and the Darcy's law (J7|). 
Thus, the complete set of equations required to fully determine the four dynamic variables h, 
p, v r and Vcf, consists of four differential equations (together with the appropriate boundary 
conditions at the contact line): one equation of the conservation of mass (|2|L two equations 
of the Darcy's law (J7j), and one equation of the mechanical equilibrium of the interface Q. 
They provide all the necessary conditions to solve the problem at least in principle. 

In practice, however, solution of these four coupled differential equations is not possible in 
the geometry of interest. At the same time, under normal drying conditions the viscous stress 
is negligible, or, equivalently, the typical velocities are much smaller than v* = o~/3r) ~ 24 m/s 
(for water). As shown in the Appendix, the four equations decouple under these conditions. 
As a result, one can employ the equilibrium result for the surface shape h(r, <fi) at any 
given moment of time, and then determine the pressure and the velocity fields for this fixed 
functional form of h. Mathematically, the original system of equations can be rewritten as: 

2H = -**, (10) 
a 

V • (/i 3 VV0 = — " d t h, (11) 
P 

v = h 2 Vip, (12) 

where Ap = p — p a tm, i> = — and po and p\ are the leading and the first-order terms 

in the expansion of pressure p = po + ep\ + • • • in a small parameter e inversely proportional 
to v* (see the Appendix for details). Note that po is independent of (r, 0), although it does 
depend on time (and this time dependence will be determined later in this paper). Therefore, 
there is a profound difference between equations (jOJ) and (fTOj) : the former is a local statement, 
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with the right-hand side depending on the coordinates of a point within the angular sector, 
while the latter is a global condition of spatial constancy of the mean curvature throughout 
the drop. Equation (fTU|) defines the equilibrium surface shape for any given value of po at any 
given moment of time, and moreover, can be solved independently of the other equations. 
Thus, the procedure for finding the solution becomes significantly simplified: first find the 
equilibrium surface shape h(r, 0) from condition (jl(Jj) and independently specify the functional 
form of the evaporation rate J(r, 0), then solve equation (fTTj) for the reduced pressure ip(r, 0), 
and finally obtain the flow field v(r, 0) according to prescription (fT2")l . The next three sections 
will be devoted to the particular steps of this procedure. 



Surface shape. There are well defined equations governing the equilibrium surface 
shape h(r, 0) and the evaporation rate J(r, 0), but there is no generic method for solving these 
equations analytically in an arbitrary geometry, in particular, in the geometry of the angular 
sector under consideration. Moreover, even if one could find these exact analytic expressions, 
a second-order differential equation of the kind of Eq. (fTTj) would not be guaranteed to 
have an analytic solution in closed form for arbitrary functions h and J. Hence, a feasible 
way to proceed analytically is to seek an approximate solution that captures the essential 
physical features and is correct at least asymptotically. Here, in the geometry of an angular 
sector, the only possible locations of singularities and divergences (which normally govern the 
properties of the solution) are at the vertex of the angle (i.e. at r = 0) and at its sides (i.e. 
at = ±a/2). Therefore, the most important physical features will be correctly reflected if 
asymptotic results (as r — > and as — > ±a/2) are found analytically. 

The boundary problem for the equilibrium surface shape of the drop consists of the 
differential equation (jl(J|) and boundary conditions at the vertex and at the sides of the 
angle: 

h(0, 0) = h(r, -a/2) = h(r, a/2) = 0. (13) 

Equation (fTUj) represents the fact that the local mean curvature is spatially uniform, but 
changes with time as the right-hand side (Ap) changes during the drying process. The 
asymptotic solution to the boundary problem (JTUJ), (fT3^) was found in our earlier paper j2]. 
The result turned out to have two qualitatively different regimes in opening angle a (acute 
and obtuse angles) and can be written as 



Mr,#) = ^. (M) 



Here R(t) = a/Ap and is the only function of time in this expression; exponent v has a 
discontinuous derivative at a = n/2 and is shown in Fig. |H1 and 

if < a < f (v = 2), 

l'(o) = { " (15) 

if | < a < n (p — tt/o). 

The constant C cannot be determined without imposing boundary conditions on h at some 
curve on the side of the drop furthest from the vertex of the angle. It is restricted by neither 
the equation nor the side boundary conditions, and thus, is not a universal feature of the 
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Figure 6: Dependence of exponent v in the power law h(r) [Eq. ()14j) ] on opening angle a, 
after Ref. 0. 

solution near the vertex of the angle. The constant C can (and does) depend on the opening 
angle a. As we showed in the earlier paper, this constant must have the following diverging 
form near a = tt/2: 

C = -^-^ + C + 0(a - tt/2) (16) 

where Co is independent of a. We will adopt this form of C (with Co set to unity) for all 
numerical estimates for obtuse opening angles. 

Two different values of v corresponding to the acute and obtuse angles give rise to the 
two qualitatively different regimes for surface shape. This difference can best be seen from 
the fact that the principal curvatures of the surface stay finite as r — > for acute angles and 
diverge as a power of r for obtuse angles. This qualitative difference can be observed in a 
simple experimental demonstration, which we provided in our earlier work We refer to 
that earlier work for further details and discussion. We only note here that the asymptotic 
r — ► at the vertex of the angle actually means r <C R (which is typically of the order of a 
few millimeters for water under normal conditions), and that V/t is indeed small for r ^ R 
and can be safely neglected with respect to unity (i.e. the free surface is nearly horizontal in 
the vicinity of the tip of the angle) as was asserted earlier. 

Evaporation rate. Since the rate-limiting step in the evaporation process is diffusion of 
saturated vapor just above the liquid-vapor interface, density of vapor n obeys the diffusion 
equation. However, diffusion rapidly attains a steady state (typically in a fraction of a 
second), and therefore the diffusion equation reduces to the Laplace equation 

V 2 n = 0. (17) 

This equation is to be solved together with the following boundary conditions: (a) along 
the surface of the drop the air is saturated with vapor and hence n at the interface is the 
constant density of saturated vapor n s , (b) far away from the drop the density approaches 
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Figure 7: Illustration of analogy between evaporation rate J for a liquid drop and electric 
field E for a conductor. Consideration of the drop (or conductor) and its reflection in the 
plane of the substrate significantly simplifies the boundary problem. 



the constant ambient vapor density n^, and (c) vapor cannot penetrate the substrate and 
hence d z n = at the substrate outside of the drop. Having found density of vapor, one can 
obtain the evaporation rate J = —DVn, where D is the diffusion coefficient. 

This boundary problem is mathematically equivalent to that of a charged conductor of 
the same geometry at constant potential if we identify n with the electrostatic potential and 
J with the electric field. Moreover, since there is no component of J normal to the substrate, 
we can further simplify the boundary problem by considering a conductor of the shape of 
our drop plus its reflection in the plane of the substrate in the full space instead of viewing 
only the semi-infinite space bounded by the substrate (Fig. |7J) . This reduces the number of 
boundary conditions to only two: (a) n = n s on the surface of the conductor, and (b) n = 
at infinity. The shape of the conductor (the drop over the angular sector and its reflection 
in the substrate) now resembles a dagger blade. So, now we have to tackle the problem of 
finding the electric field around the tip of a dagger blade at constant potential in infinite 
space. 

If one decides to account for the thickness of the blade (given by doubled h(r, <fi) of 
Eq. (|14p) accurately, it becomes apparent that there is no hope for any analytical solution in 
this complex geometry. However, taking into account that near the tip V/i is very small and 
hence the thickness of the blade itself is very small, we can approximate our thick dagger 
blade with a dagger blade of zero thickness and the same opening angle (i.e. with a flat 
angular sector). In the limit r — > the contact angle 9 scales with r as (r / ' R) v ~ l and hence 
goes to zero. Thus, only the flat blade can be considered up to the main order in r. This 
approximation would not be adequate for determining the surface shape or the flow field, but 
it is perfectly adequate for finding the evaporation rate. We will discuss possible corrections 
to this result later in this subsection. 

The problem of finding the electric field and the potential for an infinitely thin angular 
sector in three-dimensional space requires introduction of the so-called conical coordinates 
(the orthogonal coordinates of the elliptic cone) and heavily involves various special functions. 
Luckily, it was studied extensively in the past (TU1 HU E21 EE] although the results cannot 
be expressed in a closed form. An important conclusion from these studies is that the r and 
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Figure 8: Dependence of exponent /i in the power law J(r) [Eqs. (J18)) and (|27|) ] on opening 
angle a, after Refs. HUH21 ESI- 



4> dependences separate and that the electric field near the vertex of the sector scales with r 
as a power law with an exponent depending on the opening angle a: 



J oc r^~ l J(4>). 

Here 

1 RV (ft* 

(19) 



J(*)« 1 ^ 



cos 0* I (96 1 * 

and and 1^(6**, 0*) are the eigenvalue and the eigenfunction, respectively, of the eigenvalue 
problem 

-L 2 Y lt (9*,<p*)= f i(ji+l)Y ll (9*,<p*) (20) 

with Dirichlet boundary conditions on the surface of an elliptic cone (degenerating to an 
angular sector as 6* — > n). In the last relation L 2 is the angular part of the Laplacian in 
conical coordinates (r, 9*, <$*). On the surface of the sector (i.e. at 9* = 7r) the relation 
between the conical coordinate 0* and the usual polar coordinate <ft is sin0 = sin (a/2) sin</>*. 
We refer to work ^2] for further details. Here we notice only that neither fi nor Y tl {9* 1 <$*) can 
be expressed in a closed analytic form; however, the exponent \x can be computed numerically 
and is shown in Fig. |H] as a function of a. Note that this exponent is lower than similar 
exponents for corresponding angles for both a wedge (a two-dimensional corner with an 
infinite third dimension) and a circular cone. Both these shapes (wedge and cone) allow 
simple analytical solutions but none of them would be appropriate for the zero-thickness 
sector. 

Despite the unavailability of an explicit analytical expression for </(</>), its analytic prop- 
erties at cf) = and at <fi = ±a/2 are quite straightforward to infer. Indeed, J(4>) is an even 
function of 0; therefore, J'(0) = (as well as any other odd derivative on the bisector) and 

J(0) = J(O) + O(0 2 ) (21) 
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for small <ft. Obviously, J(0) is positive. On the other hand, at <p = ±a/2 the leading 
asymptotic of the evaporation rate (or the electric field) is known to be (A0) -1 / 2 with expo- 
nent — 1/2 corresponding to the edge of an infinitely thin half plane in the three-dimensional 
space ^3]. (We have introduced notation A0 = a/2 — \<p\ in the previous line.) If one were 
to correct this asymptotic in order to reflect the non-zero contact angle at the edge of the 
sector, the asymptotic form at <fi = ±a/2 would have to be written as 

J(0) « J*(A0)- A (22) 

where J* is a positive constant and 

A = 1=2°. (23) 

This result corresponds to the divergence of the electric field along the edge of a wedge of 
opening angle 29 (both the drop and its reflection contribute to the opening angle, hence a 
factor of 2) However, accounting for the non-zero 9 is a first-order correction to the 

main-order result A = 1/2. This can be seen from the expression for the contact angle: 



9 = arctan 



rp \ U — l 



R 



h' (a/2) 



j, \ u—l 



x ( r) ' (24) 



For all opening angles v > 1 (except a = tt where u—l). Thus, the correction due to the 
non-zero contact angle can indeed be neglected in the main-order results, and A should indeed 
be set to 1/2. Nevertheless, we will keep the generic notation A for this exponent in order to 
keep track of the origin of different parts of the final result and in order to account properly 
for the case a = tt in addition to the range of opening angles below tt. The numerical value 
of A will be assumed to be 1/2 in all estimates. 

For the purposes of the numerical estimates only we will employ the following simplified 
form of J(4>): 

■K<P)= (?) -0 2 ■ (25) 



This model form satisfies both asymptotics (|2*Tj) and (|22|). and it allows one to avoid the 
numerical solution of the eigenvalue problem (j20j) and thus not to repeat the elaborate treat- 
ment of works fOl EH C3 E| • This form is neither exact nor the only one satisfying the 
asymptotics, and the numerical graphs based on this form should be taken with a pinch of 
salt. However we expect it to be a good approximation to the true function J(4>), and in 
order to check this we conducted all the numerical calculations for an alternative form of 
J{4>) as well: 

m = {(^X-A 12 ■ (26) 



This alternative form also satisfies both asymptotics (|2*Tj) and f!22l) . In most cases (all but 
one) the discrepancy between the numerical results based on the two model forms did not 
exceed 20%, and this fact has convinced us that at least the orders of magnitude obtained by 
this approximation are correct. We would like to emphasize that only the numerical graphs 
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based on the choice of J{4>) are affected by these simplified forms; all the analytical results 
below do not rely on a particular form of J{<p) and use only the analytical asymptotics of the 
previous paragraph. 

Thus, we will use the following expression for the evaporation rate J: 




J(r,</>) = J [-=\ J(0), (27) 



where function J {(f)) is defined in Eq. (|19|) with asymptotics (|21j) and (J22j) . Here we broke 
down the constant pre-factor into two pieces: a distance scale yA (where A is the substrate 
area occupied by the drop) and all the rest Jo (which is of dimensionality of the evaporation 
rate). Trivially, Jq is directly proportional to the difference of the saturated and the ambient 
vapor densities {n s — n^). 

The evaporation rate does not depend on time and the same form of J applies during the 
entire drying process, as the diffusion process is steady. The same is true for the total rate 
of mass loss dM/dt since 



dM 
~dT 



/ JJl + (Vh) 2 rdrd(j) « — / J rdrdcj) oc —Jo A 
J a J A 



where the integrations are over the substrate area occupied by the drop. The constancy of 
this rate during most of the drying process was also confirmed experimentally This fact 
can be used to determine the time dependence of the length scale R of Eq. (fT4*j) (and hence 
of the pressure po of Eq. explicitly, as the mass M of a sufficiently thin drop is inversely 
proportional to the mean radius of curvature R: 

M cx (29) 

where we retained only dimensional quantities and suppressed all the numerical pre-factors 
sensitive to the details of the drop shape. From the last two equations one can conclude that 

dt \RJ pA v ; 

and remains constant during most of the drying process. Hence, 

where Ri is the initial mean radius of curvature {R4 = R{0)) and £/ is the total drying time: 

*/ « y^r- (32) 

Thus, at early drying stages {t <C £/) scale R grows linearly with time; this time dependence 
will be implicitly present in the results below. However, it is very weak at sufficiently early 
times and will be occasionally ignored (by setting R w RA when only the main-order results 
are of interest. 
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Flow field. With h and J in hand we proceed by solving Eq. for the reduced 
pressure ip. Assuming power-law divergence of ip as r — ■> and leaving only the main 
asymptotic (which effectively means that we neglect the regular term dth with respect to the 
divergent one J/p), we arrive at the following asymptotically-correct expression for ip: 



J r »-3v+i 



= - =H=TZ I (33) 



where time-dependence is implicitly present via i? and the function ■0(0) is a solution to the 
following differential equation: 

/i(0) /i 3 (0) 

Here tu 2 is a combination of the previously introduced exponents: 

UJ 2 = ( yU + i)(3 Z y_ /U _ i) (35) 

(plotted in Fig. El as a function of a). Computing v according to prescription (fT2"|) . we obtain 
the depth-averaged flow field 

v = v r r + v^> (36) 



v r (r, 0) = -{3u -n- 1)- —rz J 2 (^W (37) 



with components 

J r^"" 
P y/T~ 1 R-"+ 1 
and 

«,(r, 0) = j° r _"f - fe 2 (0)?(0)- (38) 
P VA J?-^ +1 

Thus, we need to solve Eq. (J34)) with respect to 0(0) in order to know the flow velocity. 

We were not able to find an exact analytical solution to this equation; however, we 
succeeded in finding approximate solutions on the bisector (|0| <C a/2) and near the contact 
line (A0 <C a/2), which represent the two opposite limits of the range of 0. (Again, we define 
A0 = a/2 — |0|.) Near the contact line the third term on the left-hand side of Eq. ()34|) is 
negligible with respect to the other two terms, and hence the solution is 

^(0) « f%*(0J(0dti + r(<f>) f /2 J(0^ + const (A0«a/2), (39) 



where 

p{4>)= [*h-*{C)dC. (40) 

Jo 

The asymptotic of this result is 

^(0)oc (A0^O). (41) 

(1-A 2 ) h'(a/2) 
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Figure 9: Dependence of parameter u of Eq. f)H5)) [governing solution ip(<p) of Eq. ()H4j) ] on 
opening angle a. 



This asymptotic can also be inferred directly from equation (JMI), without finding its solution. 
(The constant on the right-hand side of Eq. (}39|) turns out to be unimportant compared to 
the diverging terms.) In the opposite limit, on the bisector, the second term on the left-hand 
side of Eq. (|H4|) is negligible with respect to the other two terms, and therefore the solution 
is 

Jo ojh 3 (£) 

The asymptotic of this result is 

m oc m + \ - M^j -> o). (43) 

Again, this asymptotic can be obtained directly from equation (|H4jl . without solving it. The 
reduced pressure on the bisector ip(0) is positive, and so is its second derivative ip"(0) = 
u; 2, 0(O) — J(0)/i~ 3 (0). The latter is due to the facts that ip(4>) diverges positively at the edge 
of the drop (see Eq. (j41j0 and that "the Mexican hat" shape for the pressure as a function of 
polar angle is physically unlikely in this slow process. The value of ^(0) cannot be determined 
from the original differential equation; one needs to employ an integral condition resulting 
from the equality of the total in-flux into a sector of radius r by flow from the outer regions 
of the drop and the total out-flux from this sector by evaporation: 

Jrdrdcj). (44) 

Upon simplification this condition reduces to the following equation defining the constant 
pre-f actor ip(0): 

" 2 (u 2 h 3 ((j))iP((j)) - J(0)) d(f) = 0. (45) 



ra/2 


rr rot/2 
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Figure 10: Typical behavior of the numerical solution and the analytical asymptotics of 
function ipi^) f° r two values of opening angle. 



Obviously, ip(0) is proportional to J(0)/z~ 3 (0). Thus, approximate analytical solutions to 
Eq. (|34|) are available in the two opposite limits. 

In order to compensate for the unavailability of the exact analytical solution to Eq. (J34j) 
we also approached this problem numerically. The numerical solution to Eq. (J34j) satisfying 
conditions ()45|) for ^(O) and ip'(0) = for ^'(0) was found for the two model forms (|25|) 
and ()26|) of J(<f>) and for approximately 20 different values of the opening angle. In all 
cases perfect agreement between the numerical solution and the analytical asymptotics of 
the previous paragraph was observed. Two examples of the numerical solution together with 
the analytical asymptotics are provided in Fig. ^] for opening angles 70° and 110°. Both 
were obtained for the model form (|2*5Jl of function J(0), and the obtuse-angle graph used 
choice f!16|) for constant C. Different choice of the model form for function J(<f>) did not lead 
to any significant changes of these graphs. 

Characteristic behavior of the velocity field (jHUj) is shown in Fig. ^2 for a = 70° and 
a = 110° (again, obtained for the choice (|16|) and the model form (J25j) . but very insensitive 
to the particular form of </(</>) )• Note that despite the fact that the exponent (// — v) of the 
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Figure 11: Characteristic behavior of flow field for two values of opening angle. Each arrow 
represents the absolute value and the direction of velocity v at the point of arrow origin. 



power law in r is not a smooth function of a (Fig. HlZjL the qualitative behavior of the flow 
field does not visibly change as the opening angle increases past the right angle. 

The velocity diverges near the edge of the drop. This could have been deduced directly 
from the conservation of mass (0), where the divergent evaporation rate must be compensated 
by the divergent velocity (since the free-surface height is a regular function of coordinates 
and, moreover, vanishes near the contact line). Physically, change of volume near the edge 
becomes increasingly smaller as the contact line is approached and hence the outgoing vapor 
flux must be matched by an equally strong incoming flow of liquid. 

Streamlines. Given the velocity field in the drop, we can now compute the rate of 
deposit growth at the edge of the sector. We assume that the suspended particles are carried 
along by the flow with velocity equal to the fluid velocity. Integrating the velocity field (|37J1 - 
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Figure 12: Dependence of exponent (/i — v) in the power laws v r {r) and v^r) [Eqs. (fHTjl and 
(13%)) ] on opening angle a. 



dr i' r ,„ ^ ^i J (4>) 



-L = -(3^-^- (46) 
r# ^ v P V(0) 

we find the streamline equation, i.e. the trajectory of each particle as it moves with the fluid: 



r{(j>) = ro exp 



(3^-/i-l) / 



Q / 2 HO g 

no . 



(47) 



where we assume that is positive here and everywhere below (the generalization to the case 
of negative <fi is obvious as all functions of <fi are even). Thus, r = r when = a/2, so that 
r is the distance from the terminal endpoint of the trajectory to the vertex. In the limit (J41|) 
the integral in the exponent is 

* /2 ^{i)d£ (A0) 2 /A 

and the streamline equation reduces to 

r^r (A0-vO). (49) 

The streamlines are perpendicular to the contact line (up to the quadratic terms in A<p). 
This is in good agreement with what one would expect near the edge of the drop, since the 
azimuthal component of the fluid velocity diverges at the side contact line while the radial 
component goes to zero. In the limit (|4Hjl the integral in the exponent is 

a/2 4>(t)d£ 1 , a/2 , 
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and hence the result reads 



r ^r Q Hl^ (0-0). (51) 

Here we introduced 

= 3z/-^- 1 = 1 . , 

uj 2 -k 2 /i+l-K 2 /(3z/-/i-l) 1 ' 

and 

«• = - J(0) . (53) 

/i 3 (0)^(0) 

A few observations are in order about this limit and its exponents. First of all, k 2 is always 
positive as all the factors in Eq. (|5*3|) are. _ As we explained in the previous subsection, the 
second derivative on bisector ^"(0) = u 2, ip(0) — J(0)/i -3 (0) has to be positive as well, and 
therefore k 2 < uj 2 = (fi + l)(3z/ — fi — 1). The positiveness of exponent e follows both from 
this fact (as (3z/ — fi — 1) > for all a) and from the fact that the trajectory r(0) necessarily 
has to diverge as <fi — » (as solute comes from the outer regions of the drop). 

We cannot compute k 2 and e explicitly, since we do not know ^(0). However, we can 
gain some idea of the behavior of these indices by using approximate forms of J(4>) and ip(4>)- 
Fig. EI] demonstrates the characteristic behavior of parameter function of opening 

angle, obtained numerically on the basis of the model forms ()25|) and ()26|) for function J(0). 
Similarly, Fig. ITU shows the behavior of exponent e for the same two model forms of J(<p). 
In order to obtain these plots, equation (J34j) was solved numerically for each a, and then 
■0(0) was fixed by condition (J45)) . As can be observed in these graphs, the two model forms 
of J(4>) lead to the plots of very similar shape, but shifted by approximately 15-20% for k 2 
and by no more than 10% for e in the whole range of the opening angles. Thus, we conclude 
that Figs. and provide correct estimates for the qualitative behavior and the order of 
magnitude of parameter n 2 (a) and exponent e(a), respectively. Interestingly, the exponent e 
does not possess a sharp discontinuity of the first derivative at a = it/ 2 despite the presence 
of such discontinuity in parameter k 2 . 

Typical shape of the streamlines is shown in Fig. E3for ct = 70° and a = 110°. It was 
based on the model form (J2*B*|) for function J(0), and involved the corresponding numerical 
solutions for function i/j(4>) fFig. fTU|) employed in Eq. (j47|) . This shape is practically insensitive 
to the model form of J(0), and almost an identical copy of this graph was obtained for the 
model form (|26|) . 

The distance from a point on a streamline to the bisector scales with <p as 4>r(<j)) oc <p l ~' L 
when 0^0. Since e > 1 (Fig. IT^j) . this distance increases when (f> decreases. Thus, the 
streamlines diverge away from the bisector when 0^0, and hence they do not originate 
on the bisector. An incoming element of fluid initially located close to the bisector moves 
towards this bisector, reaches a minimum distance and then veers away towards the contact 
line. One can also arrive at the same conclusion having started from equation (|46|). Indeed, 
for small (in the limit f)43j)) the ratio of the velocity components is v r jv^ = —(3b> — fi — 
1)/(4>(uj 2 — k 2 )), or Vtj,/v r = — (f)/e. The latter ratio represents the angle between a streamline 
and a coordinate line <p = <po &t any point (0o; r (0o)) on that streamline. Since e > 1, the 
absolute value of this angle is less than \(f>\, and therefore, despite the opposite sign of this 
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Figure 13: Dependence of parameter k 2 of Eqs. and [governing exponents e of 

Eq. 7 of Eq. (|S5)L and 5 of Eq. (JMJ)] on opening angle a. The two curves correspond 

to the two model forms for function J ((f)): the solid curve is based on the choice 1)25)1 . the 
dotted curve is based on the choice (J2EJ). 
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Figure 14: Dependence of exponent e in the power law r(0) [Eq. 1)51)1] on opening angle a. 
The two curves correspond to the two model forms for function J ((f>): the solid curve is based 
on the choice f[SH|l . the dotted curve is based on the choice (J2EJ). 
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Figure 15: Typical shape of the streamlines for two values of opening angle (for the same 
values as in Figs. ITUlandlllj). 
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angle, the streamline diverges away from the bisector for small <p. This tendency can also be 
observed directly in Fig. [TBI 

Another feature apparent from Fig. is the self-similarity of all the streamlines. As 
is clear from equation 1)470. the only scaling parameter of the family of streamlines is r , 
and therefore all the streamlines can be obtained from a single streamline (say, the one with 
r = 1) by multiplying its r-coordinate by different values of r . 

Note that equation ()4*7)l does not contain Jo/p, and thus it is universally correct regardless 
of the choice of Jo- Physically, this indicates that solute particles move along the same tra- 
jectories independently of how fast evaporation occurs and hence how fast the flow is. Also, 
r(0) does not depend on the choice of the pre-factor C for the obtuse angles. Indeed, as we 
mentioned previously, i/j(Q) is proportional to J(0)/i _3 (0). Thus, i/j(4>) and ip'(4>) are propor- 
tional to 1/C 3 for obtuse angles (this could also be observed directly from equation (|S4*|i). 
Therefore, the right-hand side of the equation ()47)) is independent of C. 

These general features of r(0) are reflected in the behavior of the index k, 2 (and hence 
the exponent e and all other exponents dependent on k 2 that we will introduce later). Since 
i[>(0) is proportional to J(0)/i~ 3 (0), then n 2 is independent of J(0)/i~ 3 (0) despite the explicit 
presence of this combination in its definition. Thus, index k 2 is indeed independent of the 
evaporation intensity and the constant pre-factor of the surface shape, in good agreement 
with the general observations of the previous paragraph. On the basis of equation (j43|) 
defining ?/>(0), index k 2 can be written in the form 



demonstrating its independence from the pre-factors of each function of <p. 

Solute transfer. Three time regimes. Now, given the shape of the streamlines, 
we use our knowledge of the initial distribution of the solute, namely, that the solute has 
constant concentration c everywhere in the drop at time t — 0, and compute the time it 
takes an element of fluid (moving along a streamline) to reach the contact line at distance tq 
from the vertex having started from some point (r, <fi) on that streamline. This time can be 
found by integrating either dt = rdifi/v^ or dt = dr/v r with known or v r and the relation 
between r and <j) on the streamline (Eq. ()47j) ): 




(54) 





:/2 exp (v - fi + l)(3v - fjt - l)f ( 




(55) 




where to is a combination of system parameters with dimensionality of time: 




(56) 
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Figure 16: Qualitative sketch: mutual location of streamlines (the two lines with arrows) 
and isochrones (the four numbered lines). Solute moves along the streamlines towards the 
contact line (the bold line). Shaded area is swept by an infinitesimal element of fluid between 
the two infinitesimally-close streamlines as that element moves towards the contact line. The 
isochrones are the geometric locations, starting from which the solute reaches the contact 
line at the same time. Solute from isochrone 1 reaches the contact line first; solute from 
isochrone 4 reaches the contact line last. 



Within this time all the solute that lays on the way of this element of fluid as it moves 
toward the contact line becomes part of the deposit (highlighted area in Fig. H5j) . The mass 
dm of this deposit (accumulated on the contact line between r and r + dr ) can be found 
by integrating h(r, 0) over area dA swept by this infinitesimal volume and multiplying the 
result by the initial concentration c of the solute: 



dm = c h(r, 0) rdrdcj). 



dA 



(57) 



Employing relation (}4~Tj) once again, we obtain: 



dm 



^dro 



R 



v-l 



a/2 _ 

HQ exp 



{v + 2){3v- ii- 1) 



no 



d(. 



(5f 



Dependence dm(t) can now be found by eliminating from results (J53J) and (|58j) . Since 
we use depth- averaged velocity throughout this paper, we implicitly assume that there is no 
vertical segregation of the solute. 

Exact analytical calculation of the dependence m(t) is not possible for an arbitrary start- 
ing point (r, 0) on a streamline since no analytical expression for ip(4>) is available for arbitrary 
and since integrals in Eqs. (|55jl and (|58jl cannot be computed analytically for arbitrary 
even if ■0(0) were known. However, there are two important cases that can be tackled analyt- 
ically: (a) early times, when the initial point is close to the contact line (i.e. when A0 < a/2 
and r m r ) and only the solute between that initial point and the contact line is swept 
into the edge deposit (the starting point is on isochrone 1 of Fig. or closer to the contact 
line), and (b) intermediate times, when the initial point is close to the bisector (i.e. when 
|0| -C a/2 and r r ) and virtually all the solute between the bisector and the contact line 
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is swept into the edge deposit (the starting point is on isochrone 4 of Fig. ^] or further from 
the vertex). Situations between these two limiting cases (highlighted area in Fig. ^] demon- 
strates one of them, starting points on isochrones 2 and 3 would correspond to some other) 
can be extrapolated on the basis of continuity of the results. Since our region is indefinitely 
smaller than the drop as a whole, we may treat regions (a) and (b) assuming that a negligible 
fraction of the drop has evaporated. At some later stage that we call the late-time regime, 
an appreciable fraction of the drop has evaporated, and the fluid trajectories have reached 
back into the bulk of the drop. In this late regime our asymptotic treatments are clearly not 
adequate to describe the flow as we did not specify any details of the drop geometry far from 
the vertex. Thus we cannot treat this regime by our methods, and only the properties of 
drying process at early and intermediate stages can be found from information in hand. 

Apart from the definitions based on the trajectories, the three regimes can be equivalently 
defined in terms of time t: 

• early times: t -C t , 

• intermediate times: to <C t <C tf, 

• late times: t ~ tf. 

Here to is the characteristic time defined in Eq. ()56|) (this characteristic time depends on ro) 
and tf is the total drying time defined in Eq. ()32|) . The equivalence of the definitions in terms 
of the initial position on a trajectory and in terms of time can be seen from equation (|55|). 
At early times, A0 -C a/2, and the integral in Eq. (J55|) is much less than 1 (as the integrand 
is finite near = a/2 while the integration range shrinks when A0 — > 0). Thus, t <C t at 
early times. At intermediate times, \<p\ <C a/2, and the integral diverges as <fi — > 0. Hence, 
t ^> to at intermediate times. On the other hand, the difference between the intermediate 
and the late regimes lies in their relation to the total drying time as we explained in the 
previous paragraph. 

As is clear from the definition of the intermediate-time regime, the necessary condition 
for its existence is t <C tf, which can be reduced to (r /v^4) 3_M (^o/ R) u ~ 2 <C 1 by combining 
Eqs. fJSTH) and (|3^j). Since we always consider such r that r <C R and r <C VA, this condition 
is obeyed as long as 3 — > and v — 2 > 0. While the former condition is always true (/i < 3 
for all a), the latter condition is true only for acute opening angles [y = 2 for a < vr/2). 
Thus, the intermediate-time regime is well-defined for acute angles. For obtuse opening 
angles the situation is more complicated. Index v satisfies the opposite inequality [y < 2 
for a > vr/2), and hence for obtuse angles (ro/R) u ~ 2 ^> 1 when ro -C R. Combined with 
inequality (ro/y/A) 3 ^^ <C 1 this leads to an ambiguous result for how (r /\/A) 3 ~^(ro/ R) u ~ 2 
compares to 1 and hence how t compares to tf. This result depends on the exact relation 
between \/~A and R and on the numerical pre-factor in the definition of tf that we omitted 
everywhere (since it depends on the exact shape of the drop including the unspecified regions 
outside the sector of interest). Generically, VA oc R(t)9(t), where 9{i) is the contact angle in 
the bulk of the drop, i. e. far away from the vertex. Both R and 9 depend on time; however, the 
intermediate times are characterized by t <C tf, and hence, as can be seen from equation (|3H). 
R(t) pa Ri and 9{t) pa Q i in this regime (here 0, = 0(0)). Therefore, \[A cx R{9i, and the 
necessary condition (r / V~A) 3 ~^(r / R) u ~ 2 ^ 1 can be rewritten as (ro/y/Ay^^ 1 9\~ v . 
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The exponents on both sides of this inequality are positive for obtuse opening angles, and 
hence this condition should be expected to be satisfied for not too small initial values of the 
bulk contact angle 0j. The closer to the vertex the trajectory endpoint is, the better this 
condition is obeyed. On the other hand, the larger the opening angle is, the smaller v is, 
and hence the smaller Q\~ v is (assuming 6^ < 1). Thus, the condition of applicability of the 
intermediate-time regime is obeyed worse for larger opening angles. As we explain below, at 
exactly a = tc the intermediate-time regime is indistinguishable from the early-time regime, 
and hence should not exist. 

Results for m(t, r Q ) at early and intermediate times are presented in the next two subsec- 
tions. 



Deposit growth: early times. This regime corresponds to the growth entirely due 
to the transfer of particles originally located near the contact line. The starting point of a 
particle trajectory is characterized by r « r$ and A0 <C a/2. An element of fluid A0 away 
from the contact line reaches the contact line in time t of Eq. (|55|). The mass dm swept 
to the contact line by an element of fluid (iro-long for this time is defined by Eq. In 
the limit (|41|) the inner integral in the expressions (|55|) and (I58J) is given by equation ()48|) . 
Evaluating the outer integrals in these expressions in the limit (J41j) . expressing A0 in terms 
of time and then substituting the result into the expression for mass, we finally obtain the 
mass of the deposit as a function of time: 



dm , 



h! (a/2) 



1 + A J* 



2 

1+A 



R 



v-l 



A 



h' {a/2) 



(59) 



Here c is the constant initial concentration of the solute in the drop. Note that to also depends 
on vq. Thus, at early times the deposit grows in time as a power law 



^(t,r )oct^r^ 
dr 



(60) 



where the r§ arises from the r^ +1 pre-factor and from the ro-dependence of to- Using Eq. (|56|h 



we find 



P 



+ 



H + l) 



[1- A)(l + i/) -2/x 



(61) 



1 + A V r ' 1 + A 

and plot it in Fig. as a function of opening angle (the early-time curve). 

There are two important conclusions to be drawn from this result. One is that the power- 
law exponent of time 2/(1 + A) = 4/3 is exactly the same as in the case of a round drop 
considered by Deegan et al. pQ . This should be of no surprise since close to the side of the 
angle (as well as close to the circumference of a round drop) the contact line looks locally 
like a straight line, and the solute "does not know" about the vertex of the angle or the 
curvature of the circumference. This exponent is determined entirely by the local properties 
of an infinitesimal segment of the contact line of length dro and is independent of larger 
geometrical features of the system. 

The value 2/(1 + A) = 4/3 of the exponent of time can be obtained from a very simple 
argument, relying only on the assumptions that (a) the contact line is straight, (b) the 
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Figure 17: An illustration of the derivation of the four-thirds law for a = ir. The contact 
line is normal to the plane of the figure. Length dr is along the contact line and hence not 
shown. The flow is in the plane of the figure from left to right. 



streamlines are perpendicular to the contact line, and (c) the distribution of the solute is 
uniform. Indeed, the mass of both the water and the solute is proportional to the volume of 
an element of fluid near the contact line (Fig. IT?]) : dm oc (Al) 2 dro. All this mass should be 
evaporated from the surface of this volume element in some time t. The evaporation rate (per 
unit area) scales as J oc (AZ)~ A and therefore the rate of mass loss is JdA oc (Al)~ x+1 dr . The 
time it takes this volume to evaporate can now be found as the ratio of its mass to the rate 
of mass loss: t = dm/(JdA) oc (AZ) 1+A . Thus, (A/) oc tV(i+A) and hence dm/dr oc t 2 /( 1+A ) 
as asserted. 

The other observation is the dependence on r . Since exponent (3 is always between —1 
and 0, the singularity in r is always integrable at r = 0. Physically, this corresponds to 
the statement that the vertex of the sector does not dominate the sides and that the deposit 
accumulation at the vertex is not qualitatively different from the deposit accumulation on 
the sides. 



Deposit growth: intermediate times. The starting point of a streamline in the 
intermediate-time regime lies near the bisector and is characterized by coordinates r ^> r 
and |0| <C a/2. By the time an infinitesimal element of fluid from a vicinity of the bisector 
reaches the contact line virtually all the solute in the area between the bisector and the 
contact line will be deposited at the contact line. The analysis is similar to the previous case. 
In the limit (|43|) the inner integral in the expressions (|55jl and ()58|) is given by equation (|50j). 
Calculating the time it takes an element of fluid to reach the contact line (Eq. (|55p) and the 
mass accumulated at the contact line between r and r + dr for this time (Eq. and 
then eliminating from the two results, we arrive at the dependence of mass on time: 

dm ^ rg* 1 (/. + 1)(3^-/,-1)-k 2 ~ a 
dr 1 ' 0> ~ if -i (v - (i + l){3u - /X - 1) + 1 J 2 

> - Ai + l)(3i/ - ai - ^^(O^CO)-") 1+ t -'* +1Jt3 -- " 1J . (62 ) 



to/ 

Taking into account that to also depends on r , we finally conclude that the deposit mass 
grows as a power law 

^(t,r )cxt 5 rl (63) 
dr 
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where we introduced notations for the exponent of time 



6 = 1 + 



(i/-0 + 1)(3i/-//-1) 



(64) 



and for the exponent of r , originating from both the pre-factor 7q and the r -dependence 



An important observation is that the exponent of time stays greater than one in the 
intermediate-time regime. Thus, the rate of mass accumulation dm/dt continues to grow 
with time in this regime, and the deposit mass grows faster and faster. This result has a 
simple explanation for both the early- and the intermediate-time regimes. Since the initial 
distribution of the solute is uniform, and since the solvent evaporates, the solute concentration 
at any given volume increases with time. Thus, even though the fluid and the particles move 
along the same streamlines in practically constant velocity field (assuming that R(t) rs Ri 
at sufficiently early stages), the rate of mass accumulation also increases with time, since 
portions of solution arriving at the contact line at approximately constant rate have higher 
and higher solute concentration. Note that this mechanism and this result are in good 
agreement with a general conclusion of Deegan's works that the rate of mass accumulation 
must diverge at the end of the drying process (as t — > tf) and that all the deposit must 
accumulate at the contact line by t = tf. 

Another observation is related to the exponent of ro- Since k 2 < u 2 as we showed before, 
7 > — 1. Therefore, the mass is integrable at ro = 0, and the statement of the previous 
subsection (that the deposit accumulation at the vertex is not qualitatively different from 
the deposit accumulation on the sides) continues to hold in the intermediate-time regime as 
well. Trivially, 7 < \i. 

The exponent of r must be identically zero at any time for the opening angle of exactly 
a = ii. Indeed, at a = n the contact line is just a straight line (i.e. there is no angle at all), 
and therefore there is a full translational symmetry with respect to which point of this line 
should be called "vertex." Thus, the choice of r = is absolutely arbitrary, and there can 
be no dependence on r whatsoever. 

Indices 7 and 6 are plotted in Figs. HH1 and HH1 respectively, as functions of the opening 
angle (the intermediate-time curves). The graphs are based on the result for parameter 
k 2 , and the two intermediate-time curves on each graph correspond to the same two model 
forms (J25|) and (J26|) for function J(4>) as we used in Fig. The two model forms of J(4>) 
lead to a very small deviation for 5 (less than 5%) and to a more substantial difference for 7. 
The significant relative error in exponent 7 near the value a = n is due to the fact that this 
exponent has to be identically zero at exactly a = tt, while for the model forms of J(4>) it is 
a small, but non-zero number. Thus, the absolute error is still small, but this small absolute 
error divided by the small value of the exponent leads to a large relative difference. 

Results for the two time regimes. To facilitate the comparison of the results, we plot 
the exponents for the early- and the intermediate-time regimes in Figs. ITH1 and fEU together. 
The absolute values of each exponent are smaller in the intermediate-time regime, indicating 



of t (Eq. ©): 




(65) 
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Figure 18: Exponent of distance r in the power law dm/dr (r ) [Eqs. (foTi|) and (pjH|)] as a 
function of the opening angle for the two time regimes. The early-time curve corresponds to 
the exponent (3 of Eq. ([61)1: the intermediate-time curves correspond to the exponent 7 of 
Eq. (jnHJ). The two curves for the intermediate-time exponent correspond to the two model 
forms for function J ((f)): the solid curve is based on the choice (J2o"|) . the dotted curve is based 
on the choice (|26|). 
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Figure 19: Exponent of time t in the power law dm/dr (t) [Eqs. (JoT)|) and (jEU)] as a function 
of the opening angle for the two time regimes. The early-time curve corresponds to the 
exponent 2/(1 + A) = 4/3; the intermediate-time curves correspond to the exponent 5 of 
Eq. (J64|) . The two curves for the intermediate-time exponent correspond to the two model 
forms for function J(4>): the solid curve is based on the choice (J25)) . the dotted curve is based 
on the choice 
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that the dependence on distance and time gets weaker with time. We do not have a simple 
intuitive explanation for such behavior of these exponents. 

The intersection of the exponents near a = tt on both graphs can be attributed to a couple 
of reasons. First, it should be kept in mind that the plotted results for 7 and 5 are based on 
the relatively arbitrary choice of model forms (|25|) and (|26|) for function J(0). We suspect 
that these model equations for the reduced evaporation rate become increasingly inaccurate 
for large a. For example, J((fi) oc (cos^) -1 / 2 at exactly a = tt, which is not the same as 
either model form. Second, as was explained just after Eq. (|24|). at exactly a = tt the contact 
angle 9 is not small even for rCfi, and the correction to the exponent A due to this contact 
angle (see Eq. (J2BJ)) is comparable to the value 1/2 assumed in all numerical estimates. All 
in all, we believe that this intersection of the early- and intermediate-time exponents is an 
artifact of our formalism and should not be observed in reality, since the results for the two 
time regimes should be identical at exactly a = tt. At exactly a = tt the contact line is a 
straight line (no angle) and the trajectories are perpendicular to that straight contact line. 
Thus, there should be no differentiation between the early and intermediate times, since 
this differentiation is based on how far or how close to the bisector the initial point of the 
trajectory is located, and any perpendicular to a straight line can be called a bisector. As 
we showed above, at a = tt the exponent of r must be equal to zero at any time and the 
exponent of time must be equal to 2/ (1 + A) = 4/3 at any time. 

As we did everywhere above, we also find the numerical solution for d 2 m/dtdr (t) in 
addition to the early- and the intermediate-time analytical asymptotics. We find the time 
derivative of dm/dr instead of dm/dr itself in order to demonstrate the amount of mass 
arriving at the contact line at time t rather than the total mass accumulated by the time t. 
We employ the chain rule to obtain d 2 m/dtdro on the basis of Eq. (fo3j) for t(0) and Eq. (J5Hj) 
for dm/dr Q ((p): 



d 
dt 



' dm\ 
dr ) 



d I dm \ 
d<j) \dr J 
dt 



c r ( 
T R 



u+1 





T /i 3 (0)^(0)exp 



(/i + l)(3i/-/i- 1) 



(66) 



then use the numerical result for ip(4>) fFig. lrUJ) in order to find t(<f)) (Eq. (J55)0 and d 2 m/ dtdr (<f)) 
(Eq. ()66|)) numerically, and finally create a log-log parametric plot d 2 m/ dtdr$ vs. t, as shown 
in Fig. EUJ The two curves in Fig. |2~U1 correspond to the two values of a we used earlier (70° 
and 110°). Again, the plot is based on the model form f)25j) for function J(0), but very in- 
sensitive to the particular form of this function. This plot clearly demonstrates two different 
slopes (and hence two different time regimes) of each curve. The crossover between the two 
regimes (slopes) occurs around time t ~ to (^.e. near ln(t/to) ~ 0), and the early-time slopes 
are equal for both values of the opening angle (and equal to 2/(1 + A) — 1 = 1/3 as to be 
expected from our early time results). All these numerical results are in excellent agreement 
with our analytical predictions, and the numerical values of time exponents compare very 
well with those of Fig. HH1 (which should be corrected by —1 due to the differentiation with 
respect to time in Fig. |2"U|) . 

In a similar fashion we obtain a log-log plot for d 2 m/ ' dtdr§ as a function of tq. We fix t, 
then express tq in terms of by combining Eqs. ()55|) and (jo^|) . and finally determine ro(0) 
and d 2 m/dtdr ((j)) (Eq. ) numerically on the basis of the numerical result for ipi^) °f 
Fig. [TU| The resulting log-log parametric plot d 2 m/dtdr (r ) is shown in Fig. |^I]for the two 
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ln(t/t„) 



Figure 20: Log-log plot of the numerical solution d 2 m/dtdro(t) for two values of opening 
angle (for the same values as in Figs. ITTH ITTI and IT5 |) . 



values of the opening angle (70° and 110°). The purpose of this graph is to provide a snapshot 
of the deposit growth at any given moment of time t. For small r the accumulation of the 
solute at the contact line is in the intermediate-time regime, while for large r the growth is 
in the early-time regime. The threshold between the two regimes is defined by t = to- This 
condition can be reversed by solving Eq. (JoTlj) with respect to r . The resulting value 



defines the threshold in terms of Tq (at any moment of time t): the early regime corresponds 
to tq ^> r*, and the intermediate regime corresponds to ro C r*. As can be seen from the 
numeric plot, the regimes indeed switch at r ~ r* (i.e. near ln(r /r*) « 0). The intermediate- 
time slopes are almost equal for both graphs since the intermediate-time exponent 7 (the 
upper curves in Fig. IT5J) varies very weakly with a (7 ~ —0.135 for a = 70°, and 7 ~ —0.111 
for a = 110°). Again, the numerical results are in excellent agreement with the analytical 
asymptotics, and the numerical values of exponents compare very well with those of Fig. 

3 Discussion 

The mechanism presented here does not account for a number of additional effects that 
can modify the deposition. First, the vertical distribution of the solute was assumed to be 
homogeneous throughout the drying process, which is equivalent to assuming that vertical 
mixing is intensive. This assumption is quite important, and the results are expected to 
get modified if the true velocity profile © is used instead of the vertically-averaged velocity 
distribution. Qualitatively, the surface of the drop moves faster than the near-substrate 
layers, and therefore the particles near the surface reach the contact line sooner than those 
closer to the substrate. For example, if all solute particles are confined at the free surface, 




(67) 
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ln(r /r*) 

Figure 21: Log-log plot of the numerical solution d 2 m/dtdro(ro) for two values of opening 
angle (for the same values as in Figs. fTUl ITT| ll5landl2*U |) . Parameter r* is defined by Eq. (|67|). 



then the relevant velocity is the surface velocity. From Eq. (jUJ), this surface velocity is 3/2 
times the average velocity v(r, <fi) of Eq. (j7J). The result is to multiply the characteristic time 
to by 2/3 in the formulas above. Additional effects can be caused by gravity or convection, 
both leading to the non- uniform vertical distribution of solute. 

Second, for higher- viscosity liquids viscous stresses become more important, and the typ- 
ical flow velocities can become comparable with v* = cr/3r]. The velocity diverges at the 
edge, and condition v <^ v* does get violated at some distance from the edge. For water 
this distance is comparable to the typical size of the solute particles (which were 0.1-1 fim 
in diameter in the experiments), for higher- viscosity liquids it may become comparable to 
the size of the drop (a few millimeters). In the latter case the surface shape cannot be as- 
sumed to be equilibrium, as was assumed in this paper. Also, the characteristics of the flow 
are significantly altered if the concentration of the solute becomes large (this changes the 
"effective" viscosity of the solvent). 

Third, the size of the solute particles can introduce additional effects. Thus, if the particles 
are small, diffusion becomes important, and diffusive currents lead to the re-distribution of 
the solute. If the particles are too large, sedimentation may affect the amount of the material 
reaching the edge. Finally, we ignore the temperature effects on surface tension (Marangoni 
flows), that have also been observed in the experiments. 

Despite all these shortcomings, we believe that our theory captures the essential mech- 
anism of the deposit growth. The mechanism described is responsible for the transfer of 
the entire solute to the edge of the drop, and the effects above should not alter our main 
conclusions. 

Our results are similar in many ways to the results for the round drops in works of Deegan 
et al. [HE]. For example, the flow is capable of 100% transfer of the solute to the contact 
line, there are several distinctive time regimes, the deposit mass grows as a power law of 
time, and even the exponent of this power law is the same in the early-time regime. At the 
same time, there are a number of new features that did not exist in the round-drop case. 
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One of these features is the existence of the third time regime in addition to the two of 
the round-drop case. The intermediate-time regime for the angular sector mathematically 
corresponds to the late-time regime for the round drop, while the late-time regime for the 
angular sector does not have any analog in the round-drop case. The existence of this new 
regime is due to the fact that an angle, as a mathematical object, is infinite, while a circle 
always occupies finite area. A real drop with an angle is also a finite object and hence it 
must always have a section of the contact line connecting the two sides of the angle. This 
causes the existence of the third time regime determined by the solute coming from that far 
region and influenced by the presence of this new section of the contact line. The late-time 
results for the angular sector heavily depend on the shape of this new section. 

Probably the most exciting feature of the angular-sector solution is its dependence on 
the opening angle. Unlike the round-drop case, there is an extra free parameter of the 
problem — the opening angle of the sector. All the results, including the exponents of 
the power laws, depend on this opening angle. Note that these exponents are universal, 
i.e. they do not depend on any other parameters of the system, except for the opening 
angle. They are as universal as the exponent of distance —2 in the Coulomb's law. At the 
same time, the only parameter they depend on is extremely easy to control — preparing an 
evaporating drop one can adjust the opening angle of the contact line at his will without any 
technical elaborations. Thus, for example, by suitably choosing the opening angle (and the 
time regime), one can create a predetermined power-law distribution of the solute along the 
contact line with virtually any exponent of distance between —1 and (Fig. lTHJl . In principle, 
this feature may have significant practical applications for all the processes mentioned in the 
Introduction. Further control over the line deposition may be achieved by altering the contact 
line shape from a straight-sided angle to a curved-sided angle. 

Apart from the deposition along a line, a similar virtue of the present flow phenomenon is 
in setting up a well-defined concentration profile, also depending on a. Initially the amount 
of solute per unit area is simply proportional to the thickness h of the drop. As time goes 
by and the drop dries, the concentration profile changes, and the relative amount of solute 
per unit area increases as r decreases. This tendency gets further enhanced with time. This 
solute concentration profile may be important if one wants to deposit material within the 
drop. 

While our study does not completely explain the pointed shapes observed in Deegan's 
experiments (Fig. |2J), it does show that the deposition is strongest near the tip. Thus, the 
deposition tends to prevent de-pinning, and it is strongest where the de-pinning force is 
greatest. Hence it helps to maintain the angular shapes of Fig. |21 Our mechanism does not 
explain the particular angle observed in the experiments, but, clearly, our kind of analysis is 
a step towards understanding the chosen angles, the observed separation between the vertices 
and the deposition profile within the sector. 

Our findings have a unique signature that can be readily verified in experiments. One 
strong consequence of our theory is that the rate of increase of dm/dro has a sharp change 
of behavior as a function of r for any given time (Fig. |2"T|) . For small r (intermediate 
times) function d 2 m/dtdr (r ) varies weakly, while for large Tq (early times) it falls off more 
dramatically with increasing r . The crossover point r* (Eq. (|67j0 moves outwards as a 
power of time t with exponent — [i + 1). (Note that this exponent involves only the 
accurately known functions of a.) This crossover point and its outward motion provide a 
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clear-cut signature of our mechanism, and this signature should be the strongest for small 
opening angles (as Fig. ITH1 suggests). In order to avoid possible non-universal effects from 
late times, one needs to measure the system before the late-time regime. This measurement 
can be done by following particles in the flow, or by looking at the build-up of fluorescence 
at the contact line (both methods were used by Deegan et al. [H E3 IH E|). In order to avoid 
the uncertainties with depth-averaging, one can employ surface- confined tracer particles. 

One open question of our work is related to the unavailability of the exact form of </(</>) 
(as discussed after Eq. (|19p). Neither full analytical nor exact numerical form was available 
explicitly, and we had recourse to analytical asymptotics and approximate numerical expres- 
sions. Finding J(<f)) may be a formidable task, but can be accomplished at least in principle, 
as the earlier works on the subject suggest [101 HH E21 UHj- So, one way of determining the 
exponents more precisely is to try to determine J(0). Another way is related to creating such 
evaporating conditions that function J(0) is simpler, for instance, J(0) is just a constant. 
The latter case of the uniform evaporation rate is significantly easier to treat analytically, 
although there are questions on its experimental realization. Some further efforts may be 
devoted in this direction. 

Further work is also required in order to account for the finite width of the deposition 
region along the contact line. It is observed experimentally that the solute is spread out 
over a quite broad range near the tip and the sides of the angles. We believe that the finite 
width of the deposition region is related to the finite size of the solute particles and finite 
concentration of the solute. Near the contact line the volume fraction of particles becomes 
sufficiently high to influence (slow down) the flow that carries those particles. Further efforts 
are to be devoted to this problem in the future. 

The authors acknowledge useful conversations with Todd Dupont. This work was supported 
in part by the National Science Foundation under award number DMR-9975533 and in part by its 
MRSEC Program under award number DMR-0213745. 

Appendix 

The purpose of this section is to demonstrate that for sufficiently slow flows one can employ 
the equilibrium surface shape for finding the pressure and the velocity fields instead of having 
to solve for all the dynamical variables simultaneously. We will also quantify how slow 
"sufficiently slow flows" are. 

We start from the equation of the mechanical equilibrium of the interface Q, where 
we approximate the doubled mean curvature 2H with V 2 /i. This approximation holds true 
because the free surface of the drop over an angular sector is nearly horizontal near the vertex 
of the angle, as was shown in our earlier paper [2|, and the other terms of the functional H[h] 
are unimportant. Substitution of 



V = -vV 2 h + p atm 



(68) 



into the Darcy's law (J7J) yields 



v = v*h 2 V(V 2 h), 



(69) 
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where v* = a/3rj. Upon further substitution into the conservation of mass (J2J), we obtain 

V ■ (V/i 3 V(V 2 /i)) + - + d t h = 0, (70) 
v ' P 

which, together with Eq. 1)69)1 . constitutes the full system of equations for finding h(r,<p,t) 
and v(r, 0, t). 

Now, for water under normal conditions, rj = 1 mPa • s and a = 72 mN/m. Hence, the 
velocity scale v* is of the order of 

v* = — « 24 m/s. (71) 

Obviously, this is a huge value compared to the characteristic velocities encountered in usual 
drying process. Therefore, one can develop a systematic series expansion in small parameter 
e = v/v* (where v is some characteristic value of velocity, say, 10 /im/s): 

h = h + eh 1 + -- - + e n h n + ■■■, (72) 

v = v + evi + • ■ • + e n v n + ■ • • , (73) 

and keep only the ho and vo terms at the end in order to describe the process up to the main 
order in e = v/v*. A similar expansion can also be constructed for pressure: 

P = Po + epi + ■ ■ ■ + e n p n + ■ ■ • , (74) 

where Pi are related to hi by equation (|55|): 

p = -crV 2 h +Patm, Pi = -erV 2 /ti, etc. (75) 

Physically, condition v v * is equivalent to the statement that the viscous stress is negligible. 
Let us understand what h and v physically correspond to. 

Plugging the expansions for h and v into the system (|59j) - (|70j) . one obtains a set of terms 
for each power of e, starting from e _1 and up. Equating terms of the main order in e yields 
the following two equations 



^V(V^o) = and V • (^V(V^o)J = 0, (76) 

which both can be satisfied if and only if V 2 ho is a function of time only. Writing it as 



7 2 ; PO ~ Patm I 



vho = — — = ~wr (77) 

we immediately identify this equation with the statement of spatial constancy of the mean 
curvature of the interface, which describes the equilibrium surface shape at any given moment 
of time t (i.e. we obtained equation (JTQJ) with the desired properties of po). Thus, ho is indeed 
the equilibrium surface shape, and so is h (up to the corrections of the order of v/v*). 

Repeating the same procedure for the terms of the next order in e, we arrive at another 
two equations: 

v =^ 2 V(V 2 /n), (78) 
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vV ■ (/^V(V 2 /i!)) + - + d t h = 0, (79) 

which can be seen to be equivalent to the set of equations (fTT j) - (fT2j) upon identification 
■0 = vV 2 hi = —epi/3r]. Knowing the equilibrium surface shape ho, one can solve the second 
equation above with respect to the reduced pressure and then obtain velocity v by 
differentiating the result according to the first equation. Thus, up to the corrections of the 
order of v/v*, one can first find the equilibrium surface shape h(r,4>) at any given moment 
of time, and then determine the pressure and the flow fields for this fixed functional form of 
h, as was asserted in section Theory. 
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